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Abstract — It is shown first by Adleman that deoxyribonucleic acid 
(DNA) strands could be employed towards calculating solution to 
an instance of the NP-comp lete Ham iltonian Path Problem 
(HPP). Lipton also de monstrated that Adle man's te chniques 
could be used to solve the satisfiability (SAT) problem . In t his 
paper, it is demonstrated ho w the DNA operations presented by 
Adleman and Lipton can be us ed to develop the DNA- based 
algorithm for solving the Subset Product Problem. 
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I. INTRODUCTION 

In 1961, Feynman first offered bio-molecular computation, 
but his idea was not implemented by experiments for a few 
decades [1]. Adleman [2] in 1994 succeeded to solve an 
instance of the Hamiltonian path problem in a test tube, just 
by handling DNA strands. From [5], it was indicated that an 
optimal solution of every NP-complete or NP-hard problem is 
determined from its characteristics. DNA-based algorithms 
had been proposed to solve many computational problems and 
those consisted of the satisfiability problem [3], the maximal 
clique [6], three-vertex-coloring [7], the set-splitting problem 
[8],the dominating-set [9], the maximum cut problem [10] and 
the binary integer programming problem[ll]. One potentially 
significant area of application for DNA algorithms is the 
breaking of encryption schemes [12]-[13], From [14]-[15], 
DNA-based arithmetic algorithms are proposed. Furthermore 
from [16], DNA-based algorithms for constructing DNA 
databases are also offered. The intent of this work is to apply 
the biological operations in the Adleman-Lipton filtering 
model to develop the DNA-based algorithms for solving the 
subset product problem. 

The rest of the paper is organized as follows. Section II 
introduces DNA models of computation proposed by 
Adleman and his co-authors in details. Section III introduces 
the DNA program to solve the subset product problem from 
solution spaces of DNA strands. Conclusions are drawn in 
Section IV. 

II. DNA MODELS OF COMPUTATION 

In the last decade there have been revolutionary advances 

in the field of biomedical engineering particularly in 

recombinant DNA and RNA manipulating. Due to the 

industrialization of the biotechnology field, laboratory 



techniques for recombinant DNA and RNA manipulation are 
becoming highly standardized. Basic principles about 
recombinant DNA can be found in [20]-[24]. In this 
subsection we describe eight biological operations useful for 
solving the subset product problem. The method of 
constructing DNA solution space for the subset product 
problem is based on the proposed method in [18]-[19]. 

A (test) tube is a set of molecules of DNA (a multi-set of 
finite strings over the alphabet {A, C, G, T}). Given a tube, 
one can perform the following operations: 

1. Extract: Given a tube T and a short single strand of DNA, 
"s", produce two tubes +(T, s) and -(T, s), where +(T, s) is 
all of the molecules of DNA in T which contain the strand 
"s" as a sub-strand and -(T, s) is all of the molecules of 
DNA in T which do not contain the short strand "s". 

2. Merge: Given tubes Ti and T 2 , yield u(X 1; T 2 ), where u(Ti, 
T 2 ) = T\ U T 2 . This operation is to pour two tubes into one, 
without any change of the individual strands. 

3. Amplify: Given a tube T, the operation Amplify (T, T\, T 2 ), 
will produce two new tubes Xiand T 2 so that T\ and T 2 are a 
copy of T (Ti and T 2 are identical) and T becomes an empty 
tube. 

4. Append: Given a tube T and a short strand of DNA," s", the 
operation will append the short strand, "s ", onto the end of 
every strand in the tube T. It is denoted by Append (T, s). 

5. Append-head: Given a tube T and a short strand of DNA, 
"s", the operation will append the short strand, "s", onto 
the head of every strand in the tube T. It is denoted by 
Append-head (T, s). 

6. Detect: Given a tube T, say 'yes' if T includes at least one 
DNA molecule, and say 'no' if it contains none. It is 
denoted by Detect (T). 

7. Discard: Given a tube T, the operation will discard the tube 
T. It is denoted by Discard (T). 

8. Read: Given a tube T, the operation is used to describe a 
single molecule, contained in the tube T. Even if T contains 
many different molecules each encoding a different set of 
bases, the operation can give an explicit description of 
exactly one of them. It is denoted by Read (T). 

III. THE DNA ALGORITHMS FOR SOLVING THE 
SUBSET PRODUCT PROBLEM 



CISME Vol. 1 No. 2 2011 PP. 15-20 www. jcisme. org © World Academic Publishing 

15 



Communications in Information Science and Management Engineering 



risviK 



A. Definition of the Subset Product Problem 

Assume that a finite set S is {s 1; s 2 ,. .., s q }, where s m is the 
mth element for 1< m < q. Also suppose that every element in 
S and a constant M are a positive integer. The subset product 
problem is to determine whether there is a subset S' cz S such 
that the product of the sizes of the element in S' is exactly M? 
The subset product problem has been proved to be the NP- 
complete problem [25]-[26]. 

Consider that a finite set S = {1, 2, 3} and M=6. The total 
subsets for S are, respectively, 0, {1}, {2}, {3}, {1, 2}, {2, 3}, 
{1, 3} and {1,2,3}. According to the definition above of the 
Subset Product problem, the product of each subset is 
subsequently 0,1,2,3,2,6,4 and 6. So the solution of the subset 
product problem for S is {2, 3} and {1, 2, 3}.In other word, 
the product of element in S={2, 3}or in S={ 1,2,3} is equal to 
constant value 6. 

B. Construct the Solution Space ofDNA Strands for the Subset 
product Problem 

Assume that x q x q _i. . .x 2 Xi is a q-bit binary number, is 
encoded to represent one of the 2 q subsets of a q-element set S. 
From [18, 19], for every bit x^ representing the fcth element in 
S for 1< fc < q, two distinct 15-base value sequences are 
designed. One represents the value "0" for x k and the other 
represents the value "1" for x k . For the sake of convenience in 
our representation, assume that x^ 1 , which represents the fcth 
element in set S belongs to S', denotes the value of x^to be 1 
and Xfc°, which represents the fcth element in set S does not 
belong to S', denotes the value of x^to be 0. The following 
DNA-based algorithm is used to construct the solution space 
for 2 q possible subsets of a q-element set S in the subset 
product problem. 
Procedure Ink (T , q) 

(1) Append (T h x, 1 ). 

(2) Append (T 2 , x q °). 

(3) r = u(T 1; r 2 ). 

(4) For k- q -1 downto 1 

(4a) Amplify (T , T u T 2 ) 
(4b) Append (T b x,/) 
(4c) Append (T 2 , x, ) 
(4d) T = U(T 1) T 2 ) 
End For 
EndProcedure 



Lemma 1: The algorithm, Ink (T , q), is used to construct the 
solution space of 2 q possible subsets for a q-element set S. 

C. The Encoded Value for Every Element of Each Subset of a 
Finite Set 

For the purpose of appending the DNA strands that 
encode the size of the element s m in each subset S', we convert 
the element s m for 1< m < q representing the size of the 
element m as a 2n-bit binary number, s mi2n Sm,2n-i>- ■ ■> s m,2 s m,i- 
Suppose that s mi2n is the most significant bit and s ml is the 
least significant bit. For every bit s m _* , 1< m < q and 1 < fc < 
2n, from [18]-[19] two distinct DNA sequences are designed. 
One corresponds to the value "0" for s mi k and the other 
corresponds to the value "1" for s m ,jc. For the sake of 
convenience in our representation, we assume that s mi k 
denotes the value of s m ,j[ to be 1 and s m ,fc° defines the value of 
s m> j< to be 0. The following algorithm is employed to construct 
the binary value of the elements in the 2 q subsets of a q- 
element set S. 

Procedure EncodedValue (T , q, 2 n) 
(l)For m-ltoq 

(la)r 1 = + (r 0) x^)andr 2 =- 

(lb) For k- 2n downto n+1 

(lc) Append (Ti,s m>k ) 

End For 
(Id) For fc = n downto 1 

(le) Append {Ti,s m ^) 

End For 
(If) Append (T^s^+i 1 ) 
(lg) For k=2n downto n+2 
(lh) Append (T 2 , s m /) 

End For 
(li) For fc = n downto 1 

(lj) Append (T 2 ,s m _ k °) 

End For 
(Ik) T = u(Ti, T 2 ) 
EndFor 
EndProcedure 

Lemma 2: The binary value of the size for each element in 2 q 
subsets of a q-element set S can be constructed from the 
algorithm, EncodedValue(T , q, ri). 



(ToX) 



The result generated by Ink(T , q) is shown as follows: 



Tube 


The result generated by Ink (T , q) 


To 


{X 3 X 2 Xj , X 3 X 2 Xi , X 3 X 2 Xj , x 3 

x 2 x i , X 3 X 2 Xj , X 3 X 2 Xj , X 3 X 2 Xi , 

0-. 
X 3 X 2 Xi } 


Ti 





T 2 






D. The Construction of a Parallel One-bit Adder 

A one-bit adder is a Boolean function that forms the 
arithmetic sum of three inputs. It includes three inputs and two 
outputs. Two of the input bits represent augend and addend, 
respectively. The third input represents the carry from the 
previous lower significant position. The first output gives the 
value of the least significant bit of the sum for augend, addend 
and previous carry. The second output gives the output carry 
transferred into the input carry of the next one-bit adder. The 
truth table of the one-bit adder is shown in Table I. 
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Suppose that two one-bit binary numbers, « m . lk and a mk , 
represent the first input (addend) and the first output (sum) of 
a one-bit adder for 1< m < q and 1< k < 2n, respectively. A 
one-bit binary number, /3 mjk , is applied to represent the second 
input(augend) of a one-bit adder. Two one-bit binary numbers, 
Ymjc-i and y m jt , are used to represent the third input (previous 
carry) and the second output (carry) of a one-bit adder 
respectively. From [18]-[19], two distinct DNA sequences are 
designed to represent the value "0" and "1" for every 
corresponding bit. For the sake of convenience in our 
representations, assume that p m ^ contains the value of p m ^ to 
be 1 and fi m ^ contains the value of /? m ^ to be 0. Also suppose 
that a m -i^ denotes the value of « m _i_k to be 1 and a m .\x 
defines the value of o^.i^ to be 0. Similarly, suppose that a m ^ 
contains the value of a m ^ to be 1 and a m ^° denotes the value 
of « mk to be 0. y m k-i denotes the value of y m ^ _i to be 1 and 
Tm,k-i° contains the value of Xm>i,to be 0. y m ^ defines the value 
of YmM to be 1 and y m ^ contains the value of y m ^ to be 0. The 
following algorithm is proposed to perform the Boolean 
function of a parallel one-bit adder. 

TABLE I 
THE TRUTH TABLE OF THE ONE-BIT ADDER 



Augend 


Addend 


Previous 


Sum 


Carry 


bit 


bit 


Carry bit 


bit 


bit 


U 




















1 


1 








1 





1 
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1 
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1 





1 
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1 
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1 


1 


1 


1 


1 



Procedure ParallelOneBitAdder(T , # m -i,)( . Pm,k> 7m,k-i> m > k) 

(1) Ti= + (T , a m .i t k) and T 2 = -( T , a^/) 

(2) T 3 = + (T u A,/) and T 4 = - (T ly pj) 

(3) T 5 = + (T 2 , /?„,/) and T 6 =-(T 2 , A/) 

(4) T 7 = + (T 3 , y m ^) and T 8 = - (T 3 , y m ^) 

(5) T 9 = + (T 4 , y^i 1 ) and T w = - (T 4 , y^k-i 1 ) 

(6) Tn= + (T 5 , Xm.fc-1 1 ) and T i2 = - (T 5 , y m , k -i) 

(7) Ti 3 = + (T 6 , y m ,k-\) and T i4 = - (T 6 , y m , k -i) 

(8) If(Detect(T 7 ) = ="yes")then 

Append-head (T 7 , a m ^) and Append-head (T 7 , y mj k) 
Endlf 

(9) If(Detect(T 8 ) = ="yes")then 

Append-head (T a ,a m ^°) and Append-head (T 8 , y m ,k) 
Endlf 

(10) If (Detect (T 9 ) = ="yes") then 

Append-head (T 9 ,« m /) and Append-head (T 9 , y m /) 
Endlf 

(11) If (Detect (Tio) = ="yes") then 

Append-head (T 10 , a m /) and Append-head (T^y^ ) 
Endlf 

(12) If (Detect (T u ) = ="yes") then 

Append-head (T n , a m> k°) and Append-head (T n , y m ^) 



Endlf 

(13) If (Detect (T 12 ) = ="yes") then 

Append-head (T 12 , a^k 1 ) and Append-head (Ti 2 ,y m>k ) 
Endlf 

(14) If (Detect (T 13 ) = ="yes") then 

Append-head (T i3 , a m ^) and Append-head (T i3 , y m>k °) 
Endlf 

(15) If (Detect (T u ) = ="yes") then 

Append-head (T 14 , a m> u°) and Append-head (Xi 4 , y m> k°) 
Endlf 

(16) T = u(T 7 , T 8 , T 9 , r 10 , T n , T 12 , T 13 , T 14 ) 
EndProcedure 

Lemma 3: The algorithm, ParallelOneBitAdder (T , a m .y , 
Pm,k, Ym,k-i, m, k), can be applied to perform the Boolean 
function of a parallel one-bit adder. 

E. The Construction of a Parallel N-bit Adder Solution 

The parallel one-bit adder introduced in Section D figures 
out the arithmetic sum of two bits and a previous carry. 
Similarly, A binary parallel 2n-bit adder is also directly to 
perform the arithmetic sum for the two input operands of 2n- 
bit and the input carry by means of performing this one-bit 
adder 2n times. The following algorithm is proposed to 
perform the arithmetic sum for a parallel 2n-bit adder. 

Procedure ParallelAdder (T , a, (3, y, q, 2ri) 

(1) For k-2n downto 1 

(la) Append (T , « ,k°) 
EndFor 

(2) For m =1 to q 

(2a) Append (T , y m /) 
(2b)Fork=lto2n 

(2c) ParallelOneBitAdder (T , a m . hk , p mX y mM , 
m,k) 
EndFor 
EndFor 
EndProcedure 

Lemma 4: The algorithm, ParallelAdder(T , a, /?, y q, 2n), 
can be applied to perform the Boolean function to a binary 
parallel adder of 2n bits. 

F. The Construction of Right Shifter on DNA Computing 

The right shifter is applied to shift A register to the right, 
where A is a unsigned integers of 2n bits. A symbol "»" is 
used to represent the operation of right shift. Here the right 
shift must leave the sign unchanged. Assume that A can be 
represented as A^ 2 n, ^-d,2n-\, ■■■, Au for < d < 1, where the 
value for A^ k for 1 < k < 2n is 1 or 0. Suppose that the 
original value for A is represented as A 0j 2n, A 0i2n -i, ■■■, Ao, i, 
the intermediate value of the first right shift for A is 
represented as A lt 2 n, A\,2n-i, ■■■, M, i. In the processing of 
performing shift right, the bit A ,i is shifted out and is 
discarded, the bit Ai_2n is filled with a bit "0".The content of 
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Aoijs shifted into the position of A ok . lt where 2 < k < 2n. 
truth table for "»" (a right shifter) is shown as follows: 



The 



Input 


Output 


A-d-l.k 


Ad, k - 1 








1 


1 



For the sake of convenience, A^ k f or < d < 1 and 1 < k < 
2n is that the value of A d _ k is 1 and A d _ k ° for < d < land 1 < k 
< 2n is that the value of A d fc is 0. The following algorithm is 
proposed to perform the parallel right shifter. 
Procedure ParallelRightShifter(T A) 

(l)Append-head(T 1 ,Ao,2„ 1 ). 

(2)Append-head(T 2 ,A 0j2 „ ). 

(3) T = u(T b T 2 ). 

(4) For k = 2n - 1 to 1 

(4a) Amplify(r , T u T 2 ). 
(4b) Append^, A^ 1 ). 
(4c)Append(T 2 ,A 0j(< ). 
(4d) T = u(T u T 2 ). 
EndFor 

(5) For d = 1 to 1 
(5a)Append(r 0) A d , 2 „ ). 
(6) For k- 2n downto 2 

(6a) T 3 = +(T , A d -!, k 1 ) and T 4 = -(T Q , A d - h k 1 ). 
(6b) If (Detect(T 3 ) = = "yes") Then 

(6c)Append(r 3 ,A djk _ 1 1 ). 

Endlf 
(6d) If (Detect(T 4 ) = = "yes") Then 

(6e)Append(r 4 ,A djk _ 1 ). 

Endlf 
(6f) To = u(T 3 , T 4 ). 

EndFor 
EndFor 
EndProcedure 

Lemma 5: The algorithm, ParallelRightShifter(X A), can be 

employed to carry out the parallel right shifter. 

G. The Construction of a parallel binary multiplier 

Multiplication of two fixed-point binary number in 
unsigned representation is done by successive additions and 
shifting. This process is illustrated with the following 
numerical example. Let us multiply the two binary numbers 
0010 and 0011: 

2 0010 Multiplicand 

X 



3 


OOll 
0010 
0010 
0000 
0000 


Mulitiplier 


6 


0000110 


Product 



When the above process is implemented in a DNA 
computing environment, it is better to change the process 
slightly. Our fast algorithm combines the rightmost half of the 
product with the multiplier. The algorithm starts by assigning 
the multiplier to the right half of the product register, placing 
in the upper half. Depending on the value of the least 
significant bit of product is one or zero; we take the different 
operation steps in the procedure. The detailed iterations are 
shown as follows: 



Iteration 


Step 


Multiplicand 


Product 





Initial Values 


0010 


0000 0011 


1 


la: 

l=>Prod=Prod+ 

Mcand 


0010 


0010 0011 


2: Shift right 
Product 


0010 


0001 0001 


2 


la: 

l^>Prod=Prod+ 

Mcand 


0010 


0011 0001 


2: Shift right 
Product 


0010 


0001 1000 


3 


1: 0^>No 
Operation 


0010 


0001 1000 


2: Shift right 
Product 


0010 


0000 1100 


4 


1: 0^>No 
Operation 


0010 


0000 1100 


2: Shift right 
Product 


0010 


0000 0110 



In fact, if we ignore the signed bits, the length of an n-bit 
multiplicand and m-bits multiplier is a product that is n + m 
bits long. A binary multiplier is a Boolean function that 
performs the arithmetic multiplication for the multiplicand 
and multiplier with the same length n bits here. Therefore the 
variable length of product is 2n bits. The following algorithm 
is proposed to perform the Boolean function of a parallel 
binary multiplier. 

Procedure ParallelMultiplier(T Mcand, Mlier, Prod, q, 2n) 
^initialization of product 

(1) For k- 2n downto n+1 

Append (T , Prodo./t ). 
End for 

(2) For k = n to 1 

Append (T , Mlier _ k ). 
End for 
(3) For k = 1 to n 

T 3 = +(T , Prodo.i^and T 4 =-(T , Prod^ 1 ) 
IF (Detect(T 3 )='yes') then 

ParallelAdder(Xo, Mlier, Mcand, y, q, 2n) 

ParallelRightShifte r(T ,Prod) 
Endlf 
IF (Detect(T 4 )='yes') then 

ParallelRightShifte r(T ,Prod) 
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Endlf 
EndFor 
EndProcedure 

Lemma 6: The algorithm, ParallelMultiplier(X 0j Mcand, 
Mlier, Prod, q, 2n), can be employed to carry out the Boolean 
function of a parallel binary multiplier. 

H. Encode Any Given Positive Integer as a binary number in 
the Subset Product Problem 

Any given positive integer, M, can be converted as 2n one- 
bit binary numbers, M2„M 2n -i...M 2 Mi. The main advantage is 
that it is feasible for bit operations of the DNA algorithm in a 
parallel XOR operation. Cons tant(T M , 2n) is proposed to 
construct DNA strands for encoding M. The procedure is 
shown as follows: 
Procedure Constant (T M , 2n) 
(1) For k-2n downto 1 
(la) Append (T M , M k ) 

EndFor 
EndProcedure 

Lemma 7 : Any given positive integer M can be constructed 
with DNA strands from the algorithm, Constant (T M , 2n). 

I. Constructing the Parallel One-bit XOR Operation on Bio- 
molecular Computing 

The Exclusive-OR (XOR) operation of a bit for two 
Boolean variables A and B generates an output of 1 if both A 
and B have different values and if they are equal. The © 
symbol represents the XOR operation. The four possible 
combinations for the XOR operation of a bit with variables A 
and B are © = 0, © 1 = 1, 1 © = 1, 1 © 1 = 0. A truth 
table is usually used with logic operation to represent all 
possible combinations of inputs and the corresponding outputs. 
The truth table for the XOR operation is shown as follows: 



(2) T 3 = +(2\, B,, k l ) and T 4 = -(T lt B q , k ') 

(3) T 5 = +(T 2 , B„, k : ) and T 6 = -{T 2 , B q , k ') 

(4) If (Detect(T 3 ) = = "yes") then 
(4a) Append-head(T 3 , C,_ °) 
Endlf 

(5) If (Detect(T 4 ) = = "yes") then 
(5a) Append-head (T 4 , Cq^ 1 ) 
Endlf 

(6) If (Detect(T 5 ) = = "yes") then 
(6a) Append-head(X 5 , C q _ k ) 
Endlf 

(7) If (Detect(r s ) = = "yes") then 
(7a) Append-head(T 6 , C q>k °) 
Endlf 

(8) T = u(T 3 , T 4 , T 5 , T 6 ) 
EndProcedure 

Lemma 8: The procedure, 
ParallelOneBitXOR(T , A q> k , B q> k , q, k), can be 
employed to finish the parallel XOR one-bit 
Operation. 

J. Constructing the Parallel 2N-Bit XOR Operation on Bio- 
molecular Computing 

The parallel XOR operation of 2n bits simultaneously 
generates the corresponding 2n-bit outputs for XOR operation 
with two 2n-bit Boolean variables A, represented by A qt2n A qi2n . 
i..A q ,2 A q,u and B , represented by B q: 2nB q ,2n-i...B q , 2 B qil , The 
following procedure is offered to perform the parallel 2n-bit 
XOR operation for 2 q subsets of a q-element set S. 
Procedure ParallelXOR(T , A, B, q, 2n) 
(1) For k = 1 to 2n 

(la) ParallelOneBitXOR(T , A,, k , B q , k , q, k). 

EndFor 
EndProcedure 



Input 




output 


A 


B 


C=A©B 














1 


1 


1 





1 


1 


1 






Assume that two one-bit binary numbers, A q> k and B q k for 
1 < k < 2n, are used to respectively represent the first input 
and the second input for the XOR operation of a bit. Also 
suppose that a one-bit binary number, C q k for 1 < k < 2n, is 
employed to represent the output for the XOR operation of a 
bit. For the sake of convenience, assume that A q / denotes 
that the value of A q> k is 1 and A q> k ° denotes that the value of A q> 
k is 0. Similarly, suppose that B qk denotes that the value of B q 
k is 1 and B q k ° denotes that the value of B q k is 0. Assume that 
C q k denotes that the value of C q> k is 1 and C q> k ° denotes that 
the value of C qk is 0. The following procedure is proposed to 
perform the parallel one-bit XOR operation. 
Procedure ParalleOneBitXOR(T , A q k , B, k , q, k) 

(i) r 1 = +(r 0) A q , fc 1 )andr 2 = -(r 0) A (; ,' k 1 ) 



Lemma 9: The procedure, ParallelXOR(T , A, B, q, 2n), can 
be used to finish the parallel XOR operation of 2n bits. 

K. A DNA Algorithm for Solving the Subset Product Problem 

The following algorithm is used to solve the subset 
product problem of a q-element set S. Notations used in the 
following algorithm are denoted in the previous sections. 
Algorithm 1: Solving the Subset Product Problem. 

(1) Init(T , q) 

(2) EncodedValue(T , q, 2n) 

(3) ParallelMultiplier (T 0> Mcand, Mlier, Prod, q, 2n) 

(4) Constant (T M , 2n) 

(5) ParallelXOR(T , Prod, M, q,2 n) 

(6) For k = 1 to 2n 

(6a) Ti = + (T , C,, ,< ) and T 2 = -(T , C q , k °) 
(6b) T = u(T u T ) 
(6c) Discard (T 2 ) 
EndFor 

(7) If (Dectect (T ) = = "yes") then 
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(7a) Read(r ) 
Endlf 
EndAlgorithm 

Theorem 1: From those steps in Algorithm 1, the subset 
product problem for 2 q subsets of a q-element set S can be 
solved. 



L. The Complexity of Algorithm 1 

Theorem 2: Suppose that a finite set S is {siS 2j ..., s q }. The 
subset product problem for S can be solved with O (q x n) 
biological operations, O (2") library strands, O (1) tubes and 
the longest library strands, 0(qx n), where 2n is the number 
of bits for representing the value of each element in S. 



IV. CONCLUSIONS 

The subset product problem proved to be the NP- 
complete problem by restriction has been solved by a number 
of different algorithms with exponential time complexity in 
conventional silicon-based computer [26]. Here the proposed 
algorithm for solving the subset product problem is based on 
basic biological operations. The number of tubes, the number 
of biological operations, the number of memory strands and 
the longest length of memory strands, respectively, are O (1), 
O (q x n), O (2") and O (q x n). 

The presented algorithm has several advantages lying in its 
massive parallelism as described below. First, these biological 
operations needed in the proposed algorithm are 
experimentally feasible in lab level [18]-[19]. Second, the 
DNA-based algorithms of polynomial time complexity are 
proposed for searching the subset product problem. Third, the 
Adleman program [18]-[19] can be applied to generate good 
DNA sequences for the constructing solution space of our 
problem. It demonstrates that the proposed algorithm has a 
lower rate of errors for hybridization. Fourth, the contribution 
of this study is that DNA-based algorithms developed herein 
can be applied to solve addition-related problems. 
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